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Abstract 



The eigenvalue equation has been found for a Hamilton function in a form independent of the 
choice of a potential. This paper proposes a modified Fedosov construction on a flat symplectic 
lO , manifold. Necessary and sufficient conditions for solutions of an eigenvalue equation to be Wigner 

functions of pure states are presented. The 1-D harmonic oscillator eigenvalue equation in the 
coordinates time and energy is solved. A perturbation theory based on the variables time and energy 
is elaborated. 



PACS numbers: 03.65. Ca 



I Introduction 



Quantum mechanics formulated in terms of linear operators acting in a Hilbert space is a potent theory 
^ , in modern physics. However, although numerous problems have been solved in the framework of this 

OO ' mathematical model, we should be aware of its areas of weakness. The first difficulty is the lack of a 

[ direct link between quantum mechanics and classical physics. The second obstacle is the quantization 

procedure, which in its original version can only be applied to Cartesian coordinates in phase spaces of 
the type M^", 

^ : 

, It seems that deformation quantization is the formulation of quantum mechanics, which overcomes 

the above mentioned obstacles. Currently there exist versions of deformation quantization adapted to 
any symplectic or even Poisson manifold (for a review see e.g. [1]). The first version of deformation 
quantization, dedicated exclusively to the symplectic phase spaces {M.^^,uj), was proposed by Moyal [5], 
who applied the ideas of Weyl [3] , Wigner [i] and Groenewold [S] . 

'^{Jj , In this paper we deal with systems with the phase space (R^,cj), which can covered with one chart. 

The coordinates q and p represent a position and a canonically conjugated momentum respectively. The 
symplectic form cu = dq A dp. Observables are assumed to be smooth real functions defined in . 

As a *-product we have chosen the Moyal product defined by the formula (compare [S]-[7]) 

A*B = Aexpf-— A,B eC°°{R'^), (1.1) 
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where the Poisson operator 'P :— ~ 'Op'dq' ^'''^ using the sign convention compatible with 

Fedosov [8], [9] so 

ih < — > 

A* B ^ A - B - —A V B + ... 

The differential definition (|l.ip of the Moyal product follows from the Fourier representation of the 
*-product [TO] 

(a * b) iq,p) = dq'dp'dq"dp"A{q',p')B{<^',p") cxp |{(g' - q){p" - p) - {q" - q){p' - p)] 

(1.2) 
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The Moyal product is closed i.e. for all functions, for which the integrals exist 

/ dqdpA{q,p) * B{q,p) = / dqdpB{q,p) * A{q,p) = / dqdpA{q,p) ■ B{q,p). 

7r2 Jr2 Jr2 

By the Moyal bracket we mean the mapping C°°{R^) x C°"{R^) C°°{R^) defined as a factor of 
noncommutativity in the Moyal product 

{A,B}m ^(a* B - B * a). (1.3) 
in\ / 

In deformation quantization the counterpart of a density matrix g is the quasi-probability measure in 
referred to as the Wigner function W. This function is related to the density matrix by the Weyl 
correspondence (see [7]). Therefore 
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(1.4) 



In the case of a pure state, when g ^ |\I')(4'|, the Wigner function equals 



or, equivalently, 



Wi^q,p) = ^ r * (. + i) * (, - i) exp (-f ) (1.5) 



1 /' + 00 , , , / ■ 



The mean value of a function A in the state represented by the Wigner function W is the integral 

(A) = / dqdpW{q,p) ■ A{q,p). (1.7) 



The eigenvalue equation for a function A is of the form 

A*Wx=XWx (1.8a) 

with the additional condition 

{A,Wx}m^O- (1.8b) 

By A we denote an eigenvalue of the function A. Since we consider only a situation in which ^ is a real 
smooth function defined on M^, its eigenvalues are real. The symbol W\ represents an eigenfunction of 
A which is also a Wigner function. Hence it is real and normalizable. In the literature you can also find 
the names: a *-genvalue equation for (|1.8a[) and a *-genfunction for Wx (see [H]). 

This paper considers the problem of finding physical solutions to Eqs. p.8a|) and (|1.8b|) . We analyze 
the case where observable A is a 1-D nonrelativistic Hamilton function with a potential V{q) such that 
the limits limg^±oo Vil) — +00. The results can be generalized in a natural way to other functions and 
the phase spaces R^", n > 1. 

There are two obstacles to be dealt with. The first one is the choice of a local chart, in which the 
formulas (|1.8ap and (|1.8bp take a 'covariant' form. We propose such a choice and present a modified 
Fcdosov algorithm to construct the eigenvalue equation directly in this chart. The 'covariant' form of the 
eigenvalue equation is the same for all Hamiltonian functions. Complete information about the potential 
V{q) is contained in symplectic connection coefficients. The new coordinates are the energy H and the 
time T. We define this time as the time of the classical motion from a turning point but there are some 
alternative choices of the coordinate canonically conjugated to energy. 

The second obstacle is a necessary and sufficient condition for a solution of (|1.8al) . ()1.8bp to be a 
Wigner function of a pure state. We propose two criteria. The first one, presented in Theorem [1] and 
Corollary [1] is based on the fact that the Wigner function of a pure state is the image of a projection 
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operator of trace 1. This condition appears in Ref. |12| . We write it in both a differential form and an 
integral form. It seems that the latter is more useful in calculations. 

Another criterion, presented in Theorem [2l for a solution of (jl.8a|) . (jl.8b|) to be a Wigner function 
of a pure state, looks similar to the necessary and sufficient condition for a Wigner function to represent 
a pure state. However, we stress that our condition can be applied to any function. This feature is 
extremely important for practical purposes. Indeed, an analysis, if a given function is a Wigner function, 
is a complicated task. This question is considered in [T2] and [13]. 

Unfortunately, we were not able to write our criteria in a covariant form. To apply them we have 
therefore to transform our solution of (|1.8a|) . (|1.8b|) into the chart {q,p). 

The example of the 1-D harmonic oscillator illustrates our considerations. One can see that in this 
case the change of coordinates {q,p) — (T, H) radically simplifies the form of the eigenvalue equation for 
the Hamiltonian. We also present a stationary perturbation theory as another example. 

We restrict our considerations to states which, in the Hilbert space formulation of quantum mechanics, 
are represented by normalizable vectors. 



II The eigenvalue equation for a Hamilton function 

In the case when the function A is Hamiltonian H ~ ^ + V{q), the system of equations (|1.8ap . (jl.Sbl) 
is of the form (compare [TS] ) 



El til, \ a V u vv E „ „ N 



r=2,4, 



where by E we denote an eigenvalue of the Hamilton function. As may be seen, we are dealing with two 
partial differential equations. An explicit form of these equations depends on the smooth potential V{q). 
The degrees of these equations are determined by V{q) and they can be infinite. 

One of methods of solving differential equations is to change the variables. We propose a special 
choice of coordinates in which the relations (|2.9ap and (j2.9bp take a covariant form. This covariant form 
highlights the geometrical nature of the eigenvalue equation. 

We will transform Eqs. (|2.9ap . ( I2.9bp into a system of equations, which locally looks the same for 
any potential. We assume that the potential V{q) is bounded from below so we put V{q) > and 
limg_>±oo y{q) = oo. It seems to be possible to extend our method to unbounded states as well. 



II. 1 The canonical coordinates: time and energy 

Our new canonical coordinates are: a time T and the energy H. We have chosen them because for high 
energies when the system behaves classically, for every state under a single constraint H ~ E ~ const, 
the probability distribution is proportional to 6{H — E). Thus, also in the quantum case, the coordinate 
H should play a dominant role. The infiucncc of the variable T canonically conjugated to energy vanishes 
in the classical limit. 

Another argument supporting this choice of coordinates has been presented in [14j . N. C. Dias and J. 
N. Prata have shown that the phase space distribution 5^,{A — a) representing the projector \a >< a\ on 
the pure eigenstate of the operator A for the eigenvalue a is a formal h- deformation of the generalized 
function 5[A — a). 

We interpret the coordinate T as a time of arrival. Its full description will be presented in the next 
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paragraph. However, formally T is a solution of the differential equation 

dT dH dTdH _ p dT dV{q) dT _ ^ 
dq dp dp dq m dq dq dp 

It is defined up to a function f{H) and it need not represent an actual time. 

A canonical transformation {q^p) — > {T,H) is nonsingular unless p = and ^ = 0. But as the 
measure of the set of these singular points equals 0, they are negligible. In our model the coordinate T 
represents the time which is necessary to reach a point {q,p) on the phase space of the system from a 
chosen turning point {qo, 0) with a fixed H ~ E ~ const. Locally 



T = m 



f 

•I On 



sgn(p) 



/go ^2mE - 2m,V{z) 
The coordinate T has an intuitive geometric interpretation. Indeed, one can see that 

dE 



sgn{p)y^2mE — 2mV{z)dz. 



(2.10) 



(2.11) 



The integral f sgn(p)^2m£; - 2mV{z)d z is an area between the phase space trajectory drawn from the 



turning point go to the point q and the position axis (see Fig. [2|). Thus time T expresses the change of 
this area relative to the change of energy. 




Figure 1: A potential V((?) function of o. 




Figure 2: The phase space diagram. 

The definition of T fails if E equals a value of potential in its local peak, because the time of arrival 
to the peak is infinite. So the time for points with reverse momentum is not defined. 
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It is possible to redefine T in such a way that the coordinate T is finite unless p = and ^ = 0. 
Moreover, as the volume of a surface of constant energy equals 0, we have no reason to modify formula 
(|2.10p . Finally, in the case when the potential V{q) has peaks, we must cover the symplectic space (M^, oj) 
with more than one chart. 



II. 2 The symplectic connection 

To construct the eigenvalue equation for Hamiltonian in the new chart (T, H) it is enough to transform 
the variables in formulas (|2.9a|) and (j2.9b| . But this operation obscures both the geometrical character of 
the change and its interpretation. Therefore, we prefer another method - writing Eqs. (|2.9a|) and (j2.9bp 
using a symplectic connection. An introduction to symplectic differential geometry can be found in |16j . 

Definition 2.1. The symplectic connection j on a symplectic manifold (W,a;) is a torsion-free 
connection satisfying the conditions 

Wy;fc = 0, 1 < i, < 2n = dimW, (2.12) 

where the semicolon ' ; ' stands for the covariant derivative. 

In the Darboux coordinates the system of equations (|2.12p reads 

^^ij:k = -llk^ij - Ijk^ii = Ink ~ lijk = (2.13) 

and ^ijk '■= "fjk^ii- Coefficients jijk are symmetric with respect to the indices {i,j,k}. The symplectic 
manifold (>V,a;) endowed with the symplectic connection 7 is called a Fedosov manifold and it is 
denoted by (W,w,7). 

Hereafter we will work in Darboux coordinates, so locally every symplectic connection 7 will be 
characterized by the coefficients jijk, which are symmetric in their indices. Local coordinates will be 
denoted by (x^, . . . , x^"). 

The general transformation rule for symplectic connection coefficients is of the form [17) 

^ , , .i.-, . |;|:|;.,„(x>. ^ , . + ^J^^. (2.4) 



Locally the symplectic curvature tensor components are defined as 

K,M^\. . . , x^") := Lo.uK^,i = ^ _ ^ + oo'^juHsjk ~ u^^'lukls.u (2.15) 
where ujkiuj^-' — Si- 

The phase space (R'^,dq A dp) is assumed to be symplectic fiat- i.e. all symplectic curvature tensor 
components (|2.15p disappear. Moreover, in the coordinates {q,p) all the coefficients of the symplectic 
connection vanish. Hence in the new coordinates (T, H) we obtain 



dq d^p dp d^q 

^''''^^'^'^WWdQ^'WWdQ^' (2.16) 



where =T,Q^ = H. 

After simple but tedious calculations we see that 

2 d'V{q) (dV{q)V l(dV{q) 

7111 = —71 — 1 — , 7112 = - — T 7111 . T 

aq^ \ dq J p \ dq \Op 



dT 



p2 ^ \dpj ^ \ dq ^^"^ 



7i22 = --|i+l-) 1^+P7ii2) I '^-2 = ^(,a^j 
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II. 3 The Fedosov construction on a flat symplectic manifold 

Originally the Fedosov algorithm was used to introduce a Weyl type ^-product on any symplectic mani- 
fold. Its complete description can be found in [8] , [9] . In this subsection we present an adaptation of the 
Fedosov construction of the ^-product to the case of a symplectically flat Fedosov manifold (yV,aj,7). 
This procedure allows computing the Moyal product in an arbitrary chart without any reference to the 
coordinates {q,p). The Einstein summation convention is used. 

Let (W, uj, 7) be a Fedosov manifold locally covered by a chart {x^, . . . , x^"). The symbols y^, . . . , y^" 
represent the components of an arbitrary vector y belonging to the tangent space TpW at the point 
p £ W with respect to the natural basis {-^)^ ■ 

We introduce the formal series of polynomials of y^, . . . , y^" 

00 {2] 

« = E E ^'«fc.n.....„ {y'r ■ ■ • iy^n''" , (2.17) 

z=0 fc=0 

where < zi, . . . , i2„ < z — 2k , ii + ■ ■ ■ + 12,1 = z — 2k , ak,ii...i2n G C. The symbol [|] denotes the 
integer part of |. The degree of an element dk.ii...i2n is the sum 2fc + ii + . . . + 12,1. The set of the formal 
series (|2.17p at the point p is denoted as Pp*yV'[[?i]]. 

In the set Pp*yV[[/i]] we define an associative o-product. Since this new multiplication is also C[[?i]]- 
bilinear, it is sufhcient to determine the values of the o-product of the elements 

(2/l)''i...(y2«y2„o(ylpi...(y2n),.„^ 

In Darboux coordinates (compare [TS]) 

min[r,fc]+min.[j,s] / t 

{yr{f+yo{yr{y'+")'' = rijisik\ ^ (^^j {yr+'-^if+n'^'-'x 

min[j, s,t] ^ 

X E i-^T-n TTTi i<n (2.18a) 

^ ^ a(t-a)(r-t + a)(i-a)(s-a)ik-t + a) ~ ^ ' 

and 

(.yroiyr = iyriyr (2.18b) 

for |i — j| 7^ n. 

The pair {P*W[[h]],o) is called the Weyl algebra. The sum V*W[[h]] := Upew(-fp ^[[ft]], o) is known 
as the Weyl algebra bundle. 

An important role in the Fedosov construction is played by differential forms with values in the Weyl 
bundle. Locally such a 1- form equals 

oo [2] 

a = E E • ■ • (2/'")'^"dx^ (2.19) 

z=0 k=0 

The elements dk,ii...i2„,si...sm (^^7 • ■ • : 2;^") are components of a smooth symmetric tensor field on W and 
C°°(TyV) 9 y '""i;"^ ^ smooth vector field on W. Henceforth, we will omit the variables in 

ji ji si s {x'^, ■ ■ ■ , x^"). The differential 1- forms (|2.19|) are smooth sections of the bundle 7'*yV[[/i]] (E) 
A\ 

The Fedosov construction requires the use of an antiderivation operator : C°° {'P*W[[h]] ^ A^) 
C°°{V*W[[h]]), which is defined as 
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where I is the degree of a in y^ 's and it equals the number of y-^ 's. The operator 6 ^ raises the degree of 
the forms from 7'*yV[[?i]] (g) in the Weyl algebra by 1. 

The exterior covariant derivative :— dx^ A of a 0- form a in a Darboux chart is of the form 

d^a ^ da -\ — (700 — 007). (2-21) 

h 

The 1-form 7 of the symplectic connection equals 7 := ■^'fijky^y^dx'^. 

Assume that the Fedosov manifold (yV,aj,7) is symplectic flat. With each oq G C°°(W) wc assign an 
element 'P*yV[[?i]] 3 a (T^^(ao) determined by the iteration 



a ^ ao + S^'^ {d^a) . (2.22) 
Hence the component a[z] of a of the degree z equals 

J a[0] = flo, 

i a[z] = d-^(dya[z-l]), z>l. 



The projection a{a) of a € 7'*>V[[/i]] on the base space W is defined as cr(a) := a|y=o = o-o- The ^-product 
of functions ao, 60 £ C°°(W) calculated according to the rule 

ao * bo ^(^-^(ao) o (J'Hbo)) ■ (2.24) 

is the Moyal product written in the chart {x^, . . . , a;-^"). 

To illustrate the modified Fedosov construction we can calculate the *- square L * L of the angular 
momentum of a particle moving on the plane XY. The Fedosov manifold of this particle is (R^, dx Apx + 
dy Apy,"f) and in the coordinates {x,y,Pa:,Py) the symplectic connection 7 disappears. By definition 
L = xpy — ypx and after long calculations performed in the chart {x,y^px,Py) we obtain the result that 
the Moyal product L * L ^ L ■ L 

We can calculate the *- square L * L directly in the new Darboux coordinates 



r = \/x'^ + y'^, 

(j) ^ arctan(l), 
Pr — / 2 , ^Px H / 2 , nPy-: 
L = xpy- ypx- 

In the new chart {r, (j),pr, L) the symplectic form equals uj = dr /\ dpr + dcj) A dL. The nonvanishing 
symplectic connection coefficients are 

2L lor 

7112 ~ -p^ ' 7122 = —Pr , 7124 = ^~ ' T222 = -2L , 7223 = r. 
Thus the symplectic connection 1- form equals 



^y^y^dcj) - -y^y'^dL - pry^y^dcj) + ^y^y'^dr + 



-^y^y^dr - Ly'^y'^dcj) + ry'^y^dcp ~ -y'^y^dcj) - -y'^y^dr. 
2 r r 



From ((2:231) we find that 



a~\L) = L - ^yiyi + prV^y^ + Ly^y^ - ry^y"^ + y^ + ^y^y\ 
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Since the product L * L is defined as the projection ayr ^{L) o a ^{L)j , we see that in fact only three 
kinds of terms contribute to the final result: L, —ry'^y^ and ^y^y'^- The product L o L = L ■ L and 

90 1 1 ^ ]2.18q[ h'^ ih ^ 'J ih n A 1 2 4 
-ryy o-yy = — T + i:ry V - irV V ^ V V V V ■ 
r ill 

But 

'^[-^ + y - " y y y y ) = 

Thus finally L*L = L- L— ^ as expected. 



Let us write the eigenvalue equation for the Hamilton function H ~ 2m ^^^'i^ ™^ coordinates 
(T, H) . The symplectic connection is calculated according to the rule (|2.16p . As the Moyal product 
is of the Weyl type, for any A, B e C°°{M?) 

00 

fc=0 

where for every k the complex conjugation Ck{A,B) = Ck{B,A) and Ck{A, B) = {—l)''Ck{B, A). So 
for odd k's ^{Ck{A,B)) = and for even fc's ^(Ck{A, B)) = 0. By Ck we denote a bilinear differential 
operator of the order k. 

In our case both the Hamilton function H and its Wigner eigenfunction We are real. Therefore the 
requirement (jl.Sbp implies 

00 

^{H *We) ^Y.^^''^^'^^k+i{H,WE) =0 (2.25) 

A:=0 

and the relation (jl.Bap reduces to 

00 

n{H *We)=Y. ^^''C2k{H, We)^E- We- (2.26) 

fe=0 



Applying the Fedosov algorithm in the chart (T, H) we find that Eqs. (|2.25p and (|2.26p turn into 

dWE ^ ^ ^ . d^+'WE „ _ . 

^+ 2: 2: ^ (2-27a) 

r=3,5,... s+t=l 



i?.W-s+ E E ^'0-*£S = ^^^- (2.27b) 



and 

r 

r=2,4,... s+t=l 

A simple but tedious analysis shows that the coefficients Qrst are polynomials in the symplectic connec- 
tion coefficients ^ijk and their partial derivatives. The degrees of these polynomials do not exceed r. In 
the coefficient Qrst partial derivatives of the symplectic connection coefficients of the degrees less than 
r — 1 appear. For example 

1 1 1 2 

©220 = ■^■7122 , ©210 = -71117222 — -(7122) • 
O 00 

Eqs. (|2.27ap and (j2.27bp look the same for any potential V{q). The information about the potential is 
completely contained in the symplectic connection. The highest present power of h and the orders of the 
equations ()2.27ap and (|2.27bp are the same as in Eqs. (|2.9ap . (|2.9bp . 
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Ill Physically acceptable solutions of an eigenvalue equation 



In this section we discuss methods of ehmination of nonphysical solutions of Eqs. (|1.8ap . (|1.8b[) . More 
information about Wigner functions and their properties can be found in [H], [T3], HHJ-ES]- 

The formulas (jl.Sp and (|1.6p have been written in the coordinates q and p. However, since W{q,p) is 
a function, we are able to deduce its properties in an arbitrary canonical atlas on the phase space R^. 
Local coordinates in charts belonging to this atlas will be denoted by Q and P. Below arc listed several 
properties of the Wigner functions of pure states. 

Property 1. A function W{Q, P) of a pure state, as a Wigner function, is real i.e. W{Q, P) = W{Q, P). 
Property 2. The integral 

I dQdPW{Q,P) = l. (3.28) 

Moreover, since W{Q, P) represents a pure state and the Moyal product is closed, the following 
equality holds. 



Property 3. 



/ 



dQdPW''{Q,P) = (3.29) 
zTrn, 



In a case where the symplcctic space (M^, oj) is covered by more than one chart, the integral dQdPf{Q, P) 
represents the number J^2 w • /. 

As is known from the Schwarz inequality, we can estimate the function W(Q, P). Indeed, 
Property 4. 

\W{Q,P)\<\. (3.30) 

This inequality implies that every Wigner function depends on h. The Dirac constant is a real positive 
parameter. From the mathematical point of view we can choose its value arbitrarily. If the function 
W{Q, P) did not depend on h, the Umit ft ^ oo of ((330)) would rcsuh in W{Q, P) = 0. 

Property 5. The Wigner function W{q,p) of a pure state is a continuous function with respect to q for 
any p and a continuous function with respect to p for any q onM.^. 

Proof 

Let us consider the function 

F(.)=/;jd^w(,+0v,(,-i), 

where 4* e L^(IR). This function is well defined for every q G K. Indeed, changing the variables in the 
integral we can write 



/ + 00 
du^{u)^{2q-u). 
-oo 



Both functions ^/(u) and '^{2q — u) are elements of i^(R) for every q so F{q) is their scalar product. 
Hence F(q) is well defined and finite for every real q. 

Moreover, 

Fiq,p) = ^ d^^ {^^)^{^- I) exp (-f 

as the Fourier transform of the function '^(q+^)'^(q — ^) is a continuous function of p for every 



q (compare [24]). Repeating the same operations for the formula (|1.6|) we conclude that the Wigner 
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function W{q,p) of a pure state is defined and finite at every point {q,p). Moreover, it is continuous with 
respect to q for every p. ■ 

Properties [TH5] are necessary but not sufficient conditions for solutions of the eigenvalue equation 
(|1.8ap . (jl.Sbp to be Wigner functions. 

Now let us consider the necessary and sufficient conditions for a function W{Q,P) to be a Wigner 
function of a pure state. 

Theorem 1. fT^ A real function W{Q,P) defined on the phase space is a Wigner function on a pure 
state if and only if 

i. dQdPW{Q,P) = 1 and 
It. WiQ,P)*WiQ,P)^^W{Q,P). 

The proof of this statement is a straightforward consequence of the fact that in the Hilbert space 
formulation of quantum mechanics pure states are represented by projection operators of trace 1. By 
definition, the projection operator is an operator which is self-adjoint and idempotent. Since the Weyl 
correspondence W^^ establishes a one-to-one relation between operators and functions in the space M^, 
we see that Theorem [T] holds. 

Unfortunately, the function W{Q, P) usually contains arbitrary great negative powers of h. Therefore 
calculation of the Moyal product W{Q, P) * W{Q^ P) using the differential formula (jl.ip can be extremely 
difficult. However, in the coordinates we can apply the integral definition (|1.2p of the ^-product. 

Hence 

Corollary 1. A necessary and sufficient condition for a real function W{Q,P) to represent a pure 
quantum state is that 



i. J^2dQdPW{Q,P) = 1 and 



— / dq'dp'dq"dp"Wiq' ,p')W{q" ,p")exp 



'^-^{{q' -q){v" -v)-{q" -q){v' -P)] 



= W{q,p). 
(3.31) 



The main disadvantage of Corollary [T] is the fact it refers to the canonical coordinates q an p. We 
can write formula (|3.3ip in an arbitrary chart but to do this we have to substitute the original variables 
q,q',q\p,p',p" for new ones: Q,Q' ,Q" , P, P' , P" . 

For every real function VF(q,p), for which the integral p.3ip exists, the imaginary part of p.31|) 
disappears. Therefore the condition (|3.31|) is equivalent to 



4- / dq'dp'dq"dp"W{q',p')Wiq",p") cos [|((g' - q){p" - p) - [q" - q){p' ^ p)] 



W{q,p). (3.32) 



The argument of the function cos can be expressed in terms of vector products. Indeed, let V 
(g,p,0),F' = (g',p',0),y" = (g",p",0),e = (0,0,1). The relation becomes 



— / u' uj"W{V')W{V") cos 



xV" + V xV' + V" xV\ 



= W{V). 



(3.33) 



The condition p.33p can be used in any chart preserving the linear structure of . 

Another version of the sufficient and necessary condition for a function defined on to be a Wigner 
function of a pure state is formulated below. It should be emphasized that it works only in the coordinates 
{q,p). It looks similar to the well known necessary and sufficient condition for a Wigner function to 
represent a pure state |19j . However, contrary to the aforementioned result, it applies to an arbitrary 
function. 
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Theorem 2. A real function W{q,p) defined on the phase spaceR'^ is a Wigner function of a pure state 
if and only if 

i. at every point {q,p) G is continuous with respect to q and with respect to p, 
a. dqdpW{q,p) = 1, 
Hi. for every qi,q2 M. there is 

£'('71, 92) = f{qi)g{q2), 

where 



Q{qi,q2) j dpW {^^^-Y^^P^ exp 



ip{qi - q2) 



For a differentiable function Q{qi, (72) 7^ the condition (j iiil ) is usually written in the form 

lng(gi,g2) ^ Q 
dqidq2 

Proof 

"=>" The function W{q,p) is a Wigner function of a pure state. Hence it is real, continuous with respect 
to q and p and it satisfies the condition (fTil ) . As is well known |19j , every Wigner function of a pure 
state fulfills (| nil ). 

"<^=" The function q can be written in another form 



9{x,y)= / dpW {x,p)exp 
Jm 



ipy 



From the assumption (Ql) we notice that for every fixed y the function W {x,p) exp is a continuous 
function of p. Therefore the integral J^^ dp W {x,p) exp is a continuous function of x (see ^^). 

Moreover, this integral is, up to a factor, the inverse Fourier transform of W {x,p) . This observation 
implies that for every x the function g{x, y) is a continuous function of y. 

The relation between coordinates {x,y) and ((71,(72) is linear. Thus g(qi,q2) is well defined at every 
point (51,92)- 

Since the function W{q,p) is real, there must be ^((71,(72) = g('?2,'Zi)- Therefore from (j iiil l 

7{qi)9{q2) = f{q2)g{qi)- 

At all points [qi, 52) at which £i((7i, 92) 7^ 0, there must be = = A 7^ 0, where A is a constant. 

In the case qi = 02, we obtain i^p-l = = A so ^ is a real number. Finally 

siqi) 9(91) 

\f{q)\^^AMq)?- (3.34) 

At all points at which £'((71,(72) = 0, we put f{qi) ~ g{q2) ~ 0. This means that if g{q,q) = 0, then the 
formula p.34p also remains true. Hence the equality p.34p holds for every (7 S M. 

Let us define a new function F(q) := — i — f(q). So 



g{qi,q2) = /Ml^l^i) • /MJ^fe) ^ - sgn(A)F(gi)F(g2). (3.35) 
From the definition of function ^((71, (72) we can sec that 

e((7o,(7o) = / dpWiqo,p) sgn(A)|F((7o)r (3.36) 

Jm 

Assume that A < 0. Then for every go G Hi the integral J^dpW{qo,p) < 0. Therefore J^2 dqdpW{q,p) < 
which is impossible according to assumption (pi]). Hence ^4 > 0. 
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As /jjrfgo J^dpW{qo,p) — 1, we see from p.36p we see that |F(go)P can be interpreted as a spatial 
probability density. 



A solution of (|3.36p with respect to F{qo) is F{qo) — y Jt^ dpW{qo,p) exp(i0). The point go can always 
be chosen in such a way that F(go) 7^ 0. 
We put 



F{q) 



cxp(i(/ 



dpW [p 



cxp 



ip{q - go) 



Fiqo) ^J^dpW{qo,p) 
The function F{q) belongs to space i^(M), because 

/ dqdpWiq,p)^ I dq\F{q)\^ = l 

so it represents a pure quantum state. ■ 

Another sufficient and necessary condition for a square integrable function to be a pure state Wigner 
function was presented in [T^]. It would also be possible to apply the Narcowich- Wigner spectrum. A 
detailed presentation of this topic can be found in fl3| . However, it should be remembered that exclusive 
analysis of the Narcowich- Wigner spectrum is insufficient to identify non-Gaussian pure states. 



IV The examples 

Probably the most spectacular example of an application of the method proposed in our paper is the 
1-D harmonic oscillator. We presented a construction of the eigenvalue equation based on a change of 
canonical coordinates for the 1-D oscillator in [17] . Here we recall briefly the main steps of this procedure 
and concentrate on a selection of physically acceptable solutions. 

The Hamilton function of the harmonic oscillator in the coordinates {q,p) is of the form 



H = 



We put the mass to = 1 and the frequency uj = 1. 



(4.37) 




Figure 3: The harmonic potential 
Applying the relation (|2.10|) to the potential illustrated on Fig. |3]we see that 

q = V2HcosT 

p = -V^HsinT, 0<H,0<T<2tt. 



(4.38) 
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This transformation is singular at the point ((7 = 0,p = 0) but this is the set of the measure 0. It is 
assumed that in the coordinates {q,p) the symplectic connection 7 disappears. 

From the transformation rule (|2.16p we find that in the coordinates (T, H) the symplectic connection 
coefficients are 

1 

7111 = -2i7 , 7112 = , 7122 = -7777 , 7222 = 0. 
Hence the symplectic connection 1-form (|2.2ip equals 

7 = -Hy^dT - ^yVdH - ^y'y^dT. 

Applying the Fedosov algorithm we find that the eigenvalue equation for the Hamiltonian is of the form 

dWp 

*ft^=0. (4.39b) 

From the second equality we notice that = 0. Thus the eigenvalue equation (|4.39ap reduces to the 
ordinary differential equation 

The function We{H) can be written in the form We{H) ~ exp (— ^) we{H). By an easy substitution 
we can see that the function we{H) is a solution of the equation 

d^WE , A 4iJ\ dwE 2h-4E 
Introducing the new variable y = ^ we can conclude that WE{y) is a solution of the differential equation 



d'^WE ^dwE h~2E ^ , , 



As may be seen in |25| . its solution is the linear combination 

^ 2M^^fh-2E h-2E \ ^ „ _n±2E ^ /'h+2E h + 2E 
WE = Ciy G — — — , — — — +C2e^?/ G 



^ 2h 2h ' " J " \ 2h 2h 

Symbols Ci, C2 denote real numbers and the function 

r( R \ 1 ^ ^ «(« + l)/3(/3 + 1) 

G(a, /3, y) := 1 + — + — ^ + . . . 

V.y 2!j/^ 

Hence the solution of (|4.40p is 

„ / 2H\f4H\-^ fh-2Eh~2E AH\ 



2H\ f AH\~^ ^ f h + 2E h + 2E 4ff\ 



We know that eigenvalues of the Hamilton function H satisfy the inequality E > 0. 
The integral dH exp {^) {^Y^ G ^) is divergent. 
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Moreover, unless —a G M, the series G{a,a,y) is divergent for any finite y. So there must he E ~ 
h{n + i), ne JV. But ^^2/"G'( ^n, —n, —y) ~ Ln{y), where the symbol Ln{y) denotes the Laguerre 



polynomial defined as L„(y) Er«=o(-l)'" n- 



Wc can see that the solutions of the eigenvalue equation of the harmonic oscillator belong to the set 
of functions 

W^„(i/)=C„exp(^-^^i„(^^^, neU. (4.43) 

At this moment, we do not know whether all values of n G A/" are physically acceptable. To answer this 
question we apply Theorem [5] 

Every function Wn{H) is real and continuous. They are also normalizable so for C„ = ■^-^f;— (compare 
[55]) the assumption (QI]) is fulfilled. We therefore concentrate on the integral 

\ 4ft j ^" \ Yn j 

As can be verified using the computer program Mathematica, this integral equals 

' expf-llUJ^Vexpf-iU.fe^ 



2"n!v^ "V 2ftr^"vvfty 2ftr^"vvfty 

It may be concluded that the function g{qi,q2) is a product of two functions depending only on the 
variables (71,(72 respectively. From Theorem [2] every function of the form (|4.43p is a physically acceptable 
eigenvalue function. This result is in agreement with the solution obtained from the Schroedinger equation 

- Y^*s((?) + y*£;((?) = E^E{q) (4.44) 

via the relation (|1.5I) . Indeed, as is well known [5S], the solutions of (|4.44p are of the form 'i/Eiq) = 
^ ^^^^^^-Z exp (^^fl) (\^) ^^'^ eigenvalues E ~ h{n + i), n = 0, 1, 



The second example is a 1st order stationary perturbation theory. Consider a given Hamiltonian 

H = Hf) + XHi{To, Hq), 

where A e K. Focusing on an area of the phase space {M?,dTQ A (IHq), in which \XHi{Tq, Hq)\ <C \Ho\i 
the new 'perturbed' variable is 

T^To + XT,iTo,Ho) + .... 

As the variables (T, H) are expected to be canonical, in the linear approximation with respect to A we 
find that the function Ti{Ho,To) must satisfy the partial differential equation 

dTijn.Ho) dHi{To,Ho) _ 
dTo dHo 

In the 1st order approximation we deal with the canonical coordinates 



T = To + ATi(i/o,To), 
H = Ho + XHi{Ho,To) 



From (|4.45p . up to linear terms in A we have 



To = T-XTi{T,H), 
Ho = H-XHi{T,H). 



(4.45) 



(4.46) 
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Further considerations will be addressed in the linear approximation with respect to A. 

In the coordinates (Tq, Hq) a flat symplectic connection was determined by the coefficients 7111, 7112, 7122, 
and 7222- During the transformation of the coordinates {Tq^Hq) — > {T,H) these connection coefficients 
change according to the rule (|2.14p . For example 



( 57111(2;, y) 



V dx 



(T,H 



Ti(r, H) 



dy 



(T,H 



Hi{T,Hy 



-37iii(T,i7) 



5Ti(x,2/) 



dx 



{T,H) 



-1112{T,H) 



dHi{x,y) 



dx 



d^H,{x,y) 



(T.H) 



dx^ 



(T,H) 



Thus the new symplectic connection 1- form is 

7'(r,i/)=7(T,i7) + A7i(r,iJ) 

and the new J- product of functions A and B reads 

AiB = A*B + X A*B. 

The second term A A*B represents a 'correction' to the undisturbed product A* B. 
We intend to solve the eigenvalue equation 

HiWE =EWe. 

To this end we represent the eigenvalue E as the series 

E = E(, + \Ei + ... 

and its eigenfunction as 

We = Weo + AM^Bi + . • . 
Inserting these series into (|4.48p we obtain the system of equations 

H ^Wsa = EqWeo, 



(4.47) 



(4.48) 



(4.49a) 



H * Wei + H*Weo = EiWeo + EqWei, 



(4.49b) 



Let us assume that the general solution of Eq. (j4.49ap has been found and its physically acceptable 
part Weo has been extracted. The relation (j4.49bp can be written as 



{H - Eo) * Wei = E^Weo - HTWeo- 
Multiplying both sides of ()4.50p by Weq* on the left-hand side we obtain 

= ^Ei Weq - Weo * {H*Weo) 
The first correction to the energy is 

El = 2TTh [ dTdHWEo{T,H){H^WEo){T,H) 

JR2 



(4.50) 



Now, since (|4.50p is a linear nonhomogeneous equation with respect to Wei and we have already found 
the general solution to the homogeneous equation (|4.49ap . we are able to derive the function Wei. 
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V Conclusions 



We have proposed a method of solving the eigenvalue equation for a 1-D Hamiltonian based on a change 
of canonical coordinates. Instead of the coordinates {q,p) we apply the coordinates time T and the 
Hamilton function H. This change of a symplectic chart leads to Eqs. (|2.27ap . (j2.27bp . which look the 
same for every Hamiltonian. Complete information concerning the potential is contained in symplectic 
connection coefhcients so the eigenvalue equation for energy is covariant under any change of potential. 

To write an explicit form of the Moyal product in the chart (T, H) wc modified the Fcdosov construc- 
tion so that it was not necessary to refer to the coordinates {q,p). 

Among all the solutions to Eqs. (|2.27a|) . (j2.27b[) there are some unphysical ones. To eliminate them we 
propose two criteria: Theorem [1] with its variant (|3.33p and Theorem [21 Unfortunately, the differential 
definition of the Moyal product is not really applicable to Wigner eigenfunctions, since they contain 
negative powers of h. 

The integral condition p.SSp works only when coordinates are linear functions of q and p. Also 
Theorem [2] has been written in the coordinates and we were unable to find its covariant form. It 

may be that difficulties in expressing the formula p.SSp and Theorem [2] in a form invariant under any 
canonical transformation are due to the fact that in quantum mechanics spatial and momenta coordinates 
are separable. 

It seems to be possible to apply the proposed method of dealing with the eigenvalue equation to other 
observables and phase spaces (R^",a;), n > 1. We are also considering the possibility of generalizing the 
Theorems [T] and [5] to non-normalizable states. 
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